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PERTURBATION SOLUTIONS FOR TRANSONIC FLOW 
ON THE BLADE -TO-BLADE SURFACE OF 
COMPRESSOR BLADE ROWS 

by Stephen S. Stahara, Denny S. Chaussee, 
and John R. Spreiter* 

Nielsen Engineering & Research, Inc. 


SUMMARY 


A preliminary investigation was conducted to establish the 
theoretical basis of perturbation techniques, with the objective 
of minimizing computational requirements associated with para- 
metric studies of transonic flows in turbomachines. The theo- 
retical analysis involved the development of perturbation methods 
for determining first-order changes in the flow solution due to 
variations of one or more geometrical or flow parameters. The 
formulation is primarily directed toward transonic flows on the 
blade-to-blade surface of a single blade-row compressor. 

Two different perturbation approaches were identified and 
studied. They are the linear perturbation equation method and 
the nonlinear differencing perturbation method, and applications 
and results of these methods for various perturbations are pre- 
sented for selected two-dimensional transonic cascade flows . 

The solutions obtained with the linear perturbation equation 
method generally indicate that good results can be anticipated 
for blade geometry perturbations, while less satisfactory results 
occur for perturbation changes in overall quantities, such as 
blade spacing and free-stream Mach number. These results illus- 
trate the primary limitation of the linear perturbation equation 
method; that is, the basic linear variation assumption fundamental 
to the technique is severely restrictive for sensitive near- 
critical, confined transonic flows. 

For flows with shock waves, it was found that a supplemental 
perturbation condition is necessary to account for shock dis- 
placement. Initial results obtained with the nonlinear differ- 
encing method for a supercritical cascade flow, accounting for 
shock displacement, display excellent agreement with corresponding 
nonlinear calculations. 


4c 

Stanford University, 


Dept, of Applied Mechanics, Stanford, CA. 



1 . INTRODUCTION 


In view of the growing capability (e.g., refs. 1-4) of 
sophisticated niomerical techniques for the calculation of tran- 
sonic flows in turbomachinery — a capability certain to improve in 
the future — it has become apparent that a need exists for methods 
capable of alleviating the limitations imposed on these codes by 
their expensive run times. This need becomes particularly obvious 
and impelling when large numbers of cases are required as when 
examining alternative configurations in a preliminary design study 
or when the aerodynamic characteristics of a completed design must 
be determined for a wide variety of operating conditions. While 
direct acceleration techniques, as those developed for relaxation 
solutions employing nonlinear sequence transformations (ref. 5) , 
grid refinement (ref. 6), or a multilevel grid method (ref. 7), 
or more efficient algorithm procedures, such as implicit rather 
than explicit methods (refs. 8-10) for time-marching solutions, 
do reduce time requirements for particular applications, they 
represent only a partial answer. What is ultimately needed are 
techniques to supplement these evolving methods which would mini- 
mize the actual n\imber of separate solutions required in a par- 
ticular application. This would be accomplished by extending, 
over some parametric range, the usefulness of each individual 
solution calculated by the more expensive base procedures. 

Consequently, the basic motivation underlying this study is 
to extend the usefulness of such numerical solutions computed for 
specific turbomachinery configurations and flow conditions with 
a view toward reducing the computational requirements now necssary. 
The nature of the present investigation is exploratory in the sense 
that selected aspects of the procedure will be investigated, 
including validity, range of application, and economy. Since 
the basic problem selected for study, as detailed in section 2 , 
encompasses a wide variety of parameters which can be perturbed, 
as well as a complete range of operating conditions, only certain 
selected perturbation problems will be examined in detail. The 
primary goal of this study is the examination and preliminary 
development of perturbation methods to determine first-order 
changes in the flow field due to variations in one or more geo- 
metrical or flow field parameters in the transonic flow on the 
blade-to-blade surface of a single row compressor rotor. The 
utlimate goal is to develop a combination of analytical pertur- 
bation procedures and associated computer codes that would 
minimize the amount and maximize the quality of predictive cal- 
culations needed for enconomic analysis and design study of 
transonic flows in turbomachinery. 
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2. ANALYSIS 


2.1 General Procedure 

The model problem chosen to test the perturbation concept 
has as its basis the formulation of Wu (ref. 11) used for approxi- 
mating the three-dimensional flow through turbomachinery stages. 
That formulation reduces to a pair of two-dimensional problems 
that are ultimately intended to be interacted to obtain a complete 
solution. The essential ideas of the decomposition and of the 
coordinate systems used are indicated in figure 1. The hub-to- 
casing analysis, indicated in the figure on the left, accounts 
for effects of blockage of the flow area by the blades as well as 
the geometry of the hub and casing, but not for any circumferential 
variations of the flow. It provides flow properties in a merid- 
ional plane as averaged from one blade to the next, as well as 
the coordinates R = R(m) of the curved axisymmetric stream 
surface upon which the blade-to-blade analysis is carried out, and 
the spacing b (m) between two arbitrarily selected adjacent 
stream surfaces. 

The problem of primary interest in this study is associated 
with the blade-to-blade surface. That problem represents the 
two-dimensional flow past a series of blades in the curvilinear 
coordinate system (m, 6) indicated on the right of figure 1 with 
the radius, R(m) , and thickness, b(m), of the stream sheet pro- 
vided by the previous solution of the hub- to-casing problem. 

Because of the more complex geometries and range of flow condi- 
tions involved, the blade-to-blade problem is the more difficult 
to solve. However, it is rich in variety of parameters which can 
be varied so that it provides an ideal test of the proposed per- 
turbation technique. Consequently, the purpose of this study is 
to develop perturbation procedures which would account for 
first-order variations in the various geometrical quantities, 
flow parameters and operating conditions for the solution on the 
blade-to-blade surface of a single stage compressor rotor operating 
at transonic speeds. Possible perturbations include: 

Geometrical 


• Maximum thickness ratio -_t 

• Thickness distribution - f (x) 

• Maximum camber ratio -_h 

• Camber distribution - g(x) 

• Blade spacing ratio - t/c 

• Blade setting or stagger angle - X 

• Leading- and trailing-edge radii - r^g, r^e 

• Channel radius distribution - R(x) 

• Channel thickness distribution - b(x) 
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Flow 


• Inlet Mach number - Mj 

• Inflow angle - 3i 

• Outflow angle - 3a 

• Rotational speed - n 

Operating Conditions 

• Pressure ratio 

• Total pressure 

• Total temperature 

Although ultimately a general procedure capable of treating 
each of these perturbations is necessary, initially we will 
examine, under various approximations, selected mambers of these 
in detail to ascertain typical behavior. Finally, we note that 
under certain flow conditions, such as choking, it is not possible 
to vary some of the parameters independently of others. This must 
be recognized and properly accounted for in determining the per- 
turbation solutions. 


2 . 2 Governing Equations 

The level of approximation chosen as the starting point for 
analysis of the model problem is the quasi- three-dimensional 
Euler equations written for the curved axisymmetric blade-to-blade 
stream surface depicted in detail in figure 2. Applied to the 
rotating curvilinear coordinate system (m,0) fixed to the blades, 
as indicated in that figure, the equations to be solved are: 
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Here, U and V represent the relative velocity components in the 
m and 0 directions. The relative total energy and relative 
total enthalpy (rothalpy) are defined by 

Ej. = E - VfiR , = H - VJ^R (2) 
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These equations agree with those in ref. (12) , and are incorporated 
in a computer program B2DATL (ref. 13) for transonic blade-to-blade 
analysis in turbomachinery. The right-hand- side terms in equations 
(1) arise from streamwise variation of radius R(m) and thickness 
b(m) of the stream sheet and are provided by a known hub-to-casing 
solution. They vanish if the equations are applied to a two- 
dimensional cascade flow. The corresponding equations for a non- 
rotating inertial system may be obtained by setting to zero. 

The boundary conditions corresponding to these differential 
equations are indicated in figure 3 for a typical single blade 
passage. These involve a uniform inlet condition at the upstream 
boundary, periodic conditions on flow quantities on the horizontal 
boundaries extending upstream and downstream from the leading and 
trailing edges of the blades, a zero normal velocity condition at 
the blade upper and lower surfaces, shock conditions satisfying 
the Rankine-Hugoniot relations at the embedded shock waves, and a 
uniform back pressure specification at the downstream boundary of 
the passage. In addition, a Kutta condition is imposed at the 
trailing edge, and conditions insuring the continuity of pressure 
and normal velocity are required along the slipstream emanating 
from the blade trailing edge. Location of the slipstream is 
initially unknown and must be found as part of the solution. 

While only a few results obtained by solving these equations 
have been reported (refs. 2, 12), and they should be considered 
as provisional in view of the relative coarseness of the finite- 
difference calculation, those developments nevertheless indicate 
a powerful emerging general capability for analyzing high-speed 
turbomachinery. Consequently, selection of this level of solution 
as the starting point of the perturbation analysis provides a 
meaningful choice in terms of furnishing a significant goal. In 
addition, it provides as well an important series of simpler 
subcases of the general problem to check the efficiency and 
relative worth of the technique at various levels of accuracy. 


2 . 3 Perturbation Concept and Methods 


The basic concept underlying the present procedure is that 
a range of solutions in the vicinity of a previously determined 
or base solution can be calculated to first-order accuracy in 
the incremental change of the varied parameter by determining a 
linearized unit perturbation solution Q defined according to 
the relation ^ 


Q 

< V 

Approximate solution for ’ 
conditions differing from 
those of the base solution 
by an amount characterized 
by A 


^^o 

^Base solution^ 
for some flow 
quantity Q 


A • {Qp} (3) 

^Linearized ~per^ 
turbation solu- 
tion for a unit 
change of A 
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Of primary importance is the fact that the relationship defined 
by equation (3) remains valid over a range of A, while the 
linearized unit perturbation solution Qp needs to be determined 
only once. This, of course, provides the basis of the effective- 
ness of the method. The significance of the quantity Qp is 
immediately recognized from the Taylor-series nature of equation 
(3) as the rate of change of the base flow solution Qq with 
respect to the particular quantity, say q, perturbed; that is, 

Qp = OQ/3q)o. 

With regard to determining Qp, two important alternative 
techniques are available. We refer to these as the linear pertur- 
bation equation method and the nonlinear differencing perturbation 
method. Although the primary emphasis of this initial study was 
on the linear perturbation equation method, we have in fact also 
employed the nonlinear differencing perturbation method. While 
all of these applications are discussed in detail in section 3, 
it is important to recognize at this point that the two methods 
differ in philosophy and that each has its own particular strengths 
and weaknesses. 


2.3.1 Linear perturbation equation met hod.- The linear perturba- 
tion equation method represents the classical approach for per- 
forming a perturbation analysis . An overview of some of the 
fluid dynamic problems analyzed with this technique is given in 
reference 14. As the name implies, the method proceeds by 
establishing and solving a linear equation for the perturbation. 
Although in the present application our interest is confined 
solely to the first-order term, the complete procedure is embedded 
in a rational approximation scheme capable of continuation to any 
order. The method proceeds by expanding the dependent variables 
in an ascending power series of the incremental change A in the 
varied parameter, inserting that representation into the full 
governing equations, and then assembling the result into a 
corresponding series of linear equations in ascending orders of 
A. Higher-order solutions in general depend on both the base 
flow plus the lower-order solutions. Determination of the appro- 
priate boundary conditions is done in a similar fashion. 

The power of the linear perturbation equation method is that 
it requires the calculation of only one nonlinear base solution. 
With that information in hand, any number of individual perturba- 
tions can then be calculated, subject to the particular governing 
linear partial differential equation and boundary conditions which 
spply* The disadvantages of this method are that each perturbation 
problem must be posed individually, including differential equation 
and boundary conditions. Furthermore, the governing differential 
equations and boundary conditions for the perturbations must 
usually be approximated in order to simplify them to a point 
where they can be solved rapidly relative to rerunning the base 
flow procedure. Moreover, the perturbation solutions themselves 


6 



may be quite sensitive to the base flow solutions which always 
enter into the perturbation problem through the differential 
equation and sometimes through the boundary conditions as well. 

The fundamental limitation of the linear perturbation equa- 
tion technique is the restriction of the range over which the 
perturbation procedure remains valid to a linear one. Since this 
characteristic depends upon the local behavior of the base flow 
with respect to the varied parameter, no general statement regard- 
ing range of validity is possible. Typical behavior for a given 
class of flows must be ascertained by checks with the base flow 
procedure. Initially unknown at the outset of an application with 
this technique, then, are the accuracy requirements imposed on the 
base solution by the perturbation procedure and the range of 
parameter variation over which the linear assumption is valid. 

2.3.2 Nonlinear differencing method .- The alternative method to 
the linear perturbation equation technique is the nonlinear 
differencing method. With this method, the perturbation solution 
per unit change of the varied parameter, Q , is determined by 
differencing two nonlinear base flow solutions removed from one 
another by some reasonable change of a particular flow or geo- 
metrical quantity. A unit perturbation solution is then obtained 
by dividing that result by the change in the perturbed quantity. 
Related solutions are determined by multiplying the unit perturba- 
tion solution by the desired parameter change and adding that 
result to the base flow solution. This simple procedure, however, 
only works directly for continuous (no shock) flows for which the 
perturbation changes do not alter the overall definition of the 
base flow. This would be the case for subcritical flow and for 
changes in free-stream conditions (Mj, etc.) where both blade 

and cascade geometry remain fixed. For those perturbations which 
change the geometry and hence the locations where flow quantities 
are determined (such as stagger angle, blade spacing, blade 
geometry) , a coordinate stretching is necessary to insure that 
the flow regions remain compatible and that the perturbation 
solution is properly defined. Similarly, for discontinuous flows, 
where shocks are present and move under the perturbed conditions , 
special coordinate straining is necessary to account for that 
movement again in order to define properly the perturbation 
solution. 

The attractiveness of the nonlinear differencing method is 
that it is not restricted to a linear variation range but rather 
replaces the nonlinear variation between two base flow solutions 
with a linear fit. Consequently, this de-emphasizes the dependence 
and sensitivity inherent in the linear perturbation equation 
method on the local rate of change of the base flow solution with 
respect to the varied quantity. For certain high-speed turbo- 
machinery applications where, as indicated in the sketch below, 
the confined transonic flow may be highly sensitive, the linear 
range of parameter variation could be sufficiently small to be of 
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solution 



Nonlinear method 


Geometric or 
flow parameter 


limited practical use. For those situations, the nonlinear dif- 
ferencing technique is more appropriate. Furthermore, other than 
the approximation of a linear fit between two nonlinear base 
solutions, the nonlinear differencing method is not restricted 
by further approximations with respect to the governing differen- 
tial equations and boundary conditions. Rather, it retains the 
full character of the original methods used to calculate the base 
flow solutions. Moreover, no perturbation differential equations 
have to be solved, only algebraic ones. 

The primary disadvantage of this method is that two base 
solutions are required for each parameter perturbation considered. 
For certain high-use applications, this requirement can be suffi- 
ciently disadvantageous to limit the utility of the method. For 
this reason, the primary emphasis of this initial study focused 
on the linear perturbation equation method. A description of the 
analysis of the linear perturbation method is provided in the 
following section. 


2.4 Linear Perturbation Equation Method of Analysis 

2.4.1 Alternate levels of approximation .- With the ultimate level 
of accuracy for solving the base flow set at the quasi three- 
dimensional Euler equations described in section 2.2, selection 
of the appropriate level at which to perform the perturbation 
analysis remains. The alternatives are: 

• Quasi three-dimensional Euler equation 

• Quasi three-dimensional full potential equation 

• Small-disturbance transonic equation 

The most obvious choice of perturbing the Euler equations is not 
a satisfactory one since the perturbation equations become more 
cumbersome than the original nonlinear equations. Furthermore, 
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although fast solution procedures are necessary to solve these 
perturbation equations, relaxation methods have not yet been fully 
developed (ref. 15) for such equations. Consequently, the per- 
turbation solutions would also have to proceed by a time-marching 
procedure, with little chance of being significantly more econom- 
ical relative to rerunning the base solution procedure again for 
the perturbed flow while using the original base flow as an 
initial guess. 

An alternate choice is to select a quasi three-dimensional 
full potential equation approach. The significant advantages here 
are that of working with a single partial differential equation 
for the dependent variable and that relaxation methods have already 
proven successful for solving such problems (refs. 16 and 17). 
Furthermore, the method is identical in accuracy to the Euler 
equation approach for shockless flows and generally a very good 
approximation for most flows in the transonic regime where shocks 
are weak; that is, pre-shock Mach numbers less than 1.3. Finally, 
for an important series of limiting cases, checks are available 
with less accurate but more fully developed theories, such as 
transonic small-disturbance theory. 

A solution procedure based on the quasi three-dimensional 
full potential approximation is not currently available. However, 
a well-documented computer code (B2DATL) based on the Euler 
equation formulation is available and presented in references 12 
and 13. Consequently, implementation of the perturbation problem 
based on the full potential approach would proceed as follows. 

A solution of the quasi three-dimensional Euler equations 
(1) and (2) , as calculated by the time-marching technique described 
in reference 12, is assumed to represent a quasi three-dimensional 
full potential solution. Required potential derivatives are then 
determined throughout the flow by appropriate differencing of the 
known velocity field. Next, linear perturbation equations and 
boundary conditions representing various geometrical and flow 
perturbations about this assumed full potential base solution are 
derived. Finally, the perturbation equation is solved numerically 
using a successive line overrelaxation (SLOR) algorithm. 

We note that included in this formulation is an important 
subset of test perturbation cases based on the small-disturbance 
transonic approximation. A well-tested and documented code 
(TSFOIL) described in reference 18 is currently available to 
generate appropriate small-disturbance base solutions. Although 
restricted in geometrical application to certain classes of two- 
dimensional cascades, this well-documented small-disturbance 
procedure can provide fast and accurate base solutions throughout 
the entire transonic range, including supercritical, choked, and 
supersonic cases useful for testing the perturbation theory. 
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2.4.2 Base flow boundary value problem-full pot ential approxi- 
mation.- In this section we describe the base flow boundary value 
problem based on the quasi -chree-dimensional full potential steady 
approximation to the Euler equations (1) and (2) . Considering 
the steady version of equation (1) , if the 9 momentum equation 
is multiplied by U, the m momentum equation multiplied by V, 
the result added, then that expression simplified by use of the 
continuity equation together with the isentropic relation 
dp = a^dp, the final result becomes 


. (u" 


a^) 



(V^- a^) 


9V 

R90 


+ VU 


^ 9U ^ 
.R90 9m j 


= d(Rb) 

Rb dm 


dm 


(4) 


Furthermore, if we assumed that in the absolute (nonrotating) 
system that a velocity potential $ for the complete flow exists 
such that 


U = $ 


m 


V = - $9 - f2R 

then the governing differential equation for $ becomes 
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a^ = a? + i (Y-1) [U^ + (V, - J^R) ^ - (^./R - JiR)^] 

I z 1 mu 

and the subscript i denotes uniform inlet conditions. 


(5) 


(6) 


(7) 


2.4.3 Disturbance potential boundary-valu _e pr oble m.- For the 
applications considered in this preliminary study, we are 
primarily concerned with nonrotating (ft=0) cascade flows. For 
those situations, a convenient alternative form to equations (6) 
and (7) is to subtract out the incoming free-stream components 
from the potential and recast without approximation the governing 
equations in terms of a disturbance potential (j> where 

= Cqj[(x cos + y sin 3^) + (|)] (8) 
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and 


X = m/C 
y = R9/C 


Here, is the free-stream velocity magnitude and the 

absolute angle on the stream surface of the uniform incoming flow. 
The governing equation for the disturbance potential <j) follows 
directly by inserting equation (8) into equations (6) and (7) 
with fi=0. The resultant equation is given explicitly as equation 
(A-10) in Appendix A. 


Appropriate boundary conditions for the differential equation 
for the disturbance potential <() are indicated in figure 4, where 
for convenience we have chosen the origin of coordinates to lie on 
a line halfway between the blade leading edges. We note that the 
primary difference between the disturbance full potential boundary 
condition formulation and the Euler equation formulation given in 
figure 3 involves conditions along the slipstream. For the dis- 
turbance potential <j) , the conditions of continuity of pressure 
and normal velocity imply a jump in potential across the slip- 
stream surface equal to the circulation around each blade. However, 
imposition of that potential jump along the actual location of the 
slipstream is unnecessary. To assure that the potential remains 
single-valued in the domain of interest, it is only necessary to 
introduce a cut along which the potential may be discontinuous and 
which extends from the blade surface to the downstream boundary. 

The location of the cut can be arbitrary and, in particular, can 
be chosen to lie along the horizontal boundaries extending from 
the trailing edges of the blades. The form of the periodic and 
cyclic conditions imposed at corresponding points on the upper 
and lower horizontal boundaries is 


Periodic: 4>m “ 

1 B 


Cyclic : 


4 >^ = 


(9) 


where the subscripts T and B denote the top and bottom bound- 
aries, respectively and A(|)t.e is ih® jump in the disturbance 
potential at the trailing edge of the blade. 


2.4.4 Perturbation disturbance potential boundary value problem . - 
The boundary value problems associated with perturbations about 
the base flow described above are established by deriving for each 
particular perturbation both the governing partial differential 
equation and associated boundary conditions. In general, both 
will depend on the varied parameter. The derivations are initiated 
on the assximption that 




( 10 ) 


11 



where (f)® is the base flow disturbance potential and the 

perturbation disturbance potential, that is, the disturbance 
potential associated with the parameter change away from its base 
value. Insertion of equation (10) into equations (6) to (8) for 
n=0 leads to the following form of the governing equation for 
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XX 
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+ E4>^ = F 

y 


( 11 ) 


where the coefficients [A,B,C,D,E,F] are in general functions of 


<P 


0 

x' 



(f)" ,(j)0 

^xy ^xx 



and are identified in Appendix A. 


2.4.5 Boundary conditions for perturbation disturbance potential . 
The boundary conditions governing the various perturbations 
associated with equation (11) are derived in a similar fashion. 

The general form of the boundary conditions for an arbitrary 
perturbation is provided in figure 5. As anticipated, the con- 
ditions for (j) ^ are more complex than the corresponding ones 
for the base flow disturbance potential (j)°. Based on the dis- 
turbance full potential formulation, the appropriate perturbation 
boundary conditions have been derived for a number of important 
parametric variations and are provided in Appendix B. 

It is of interest to demonstrate the general procedure 
followed in determining boundary conditions for a typical per- 
turbation. The basic coordinate system used is indicated in 
the sketch below. 



In this figure, F (x) , Fj^(x) denote the ordinates of the blade 
upper and lower surface, respectively, nondimensionalized by C, 
t/C is the blade spacing ratio, X the stagger angle, (x^,y^) 
are given by 
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cos X 


( 12 ) 


and (qi,3i)/ (qa/^a) the inflow and outflow velocity and 

absolute flow direction. A perturbation in the stagger angle X 
is considered. The boundary conditions for the blade surfaces, 
and the periodic and cyclic conditions are determined as follows. 

(1) Blade boundary conditions: The basic tangency con- 

dions on the upper and lower surfaces of the blade are given 
in terms of the full potential $ by the expressions 
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The corresponding relations for the disturbance full potential 
(p are 
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Considering a perturbation in stagger angle, we allow 


(14) 





A = X® 

+ AX 


that 






Xt = t 

sin 

x“ 

+ t cos 

X®AX + 

• • • = x^ + Ax 

yt = t 

cos 


- t sin 

X®AX + 

... = + Ay 


(15) 


(16) 


then. 


X + 


2 ' 2 


Xt Ax^ y^. 
^ 2 2 ' 2 


Ay. 


(17) 


13 



Expanding this result about the base positions, denoted by 
( ) ”, using a Taylor- series , we have 
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This leads to the result that 
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Introducing the following expansion 

(f) = 4)® + 4)^ 


AX + O(AX^) 


(19) 


( 20 ) 


for the potential into equation (19) and ordering the results 
according to powers of AX, we obtain the final form of the blade 
upper surface boundary conditions for the base flow disturbance 
potential 4>° and perturbation disturbance potential 4)^ 
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Similarly, for the lower blade we obtain the result 
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(2) Periodic condition: The basic periodic condition, 

equation (9) , is given by the statement 
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As before, by using the expansions in equations (16) and (20), 
Taylor-series expanding about the base flow values, and then 
ordering that result according to powers of AX, we obtain the 
following forms of the periodic conditions for the base flow and 
perturbation component 
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where the periodic condition on the base flow velocity components 


(25) 


has also been employed. 

(3) Cyclic boundary condition: The cyclic boundary con- 
dition, equation (9), which accounts for the jump in potential 
across the passage due to circulation is given by 
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In addition to the same changes associated with the periodic con- 
dition, the cyclic condition involves changes in as well 

which must be considered. Application of the same expansion 
procedure to equation (27) leads to the result that 
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Thus, the final form of the cyclic boundary conditions is 
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As mentioned previously, corresponding results for a number 
of important perturbations have been derived and are provided in 
Appendix B. 


2.4.6 Transonic sm all-disturbance formulation . - Finally , we note 
that a variety of subcases are contained within the full potential 
approximation described above for the base flow (eqs. (6) to (8)) 
and for the perturbation component (eq. (11) ) and their associated 
boundary conditions. These include two-dimensional cascade flows 
and two-dimensional small-disturbance cascade flows. In view of 
the importance of the two-dimensional small-disturbance approxi- 
mation to this initial study, we outline briefly here the simpli- 
fication of the full potential formulation to the small-distur- 
bance level . 


By considering the two-dimensional [ (db/dx) = (dR/dx) = 0] 
nonrotating = 0) form of equations (6) to (8) , and introducing 
the scaling parameters 
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appropriate to small-disturbance flows into those equations, the 
following differential equation in conservation form results for 
the transonic small-disturbance velocity potential 
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The corresponding differential equation for the perturbation 
transonic small-disturbance potential $ defined according to 




(35) 


is given in conservation form by 
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2.5 Base Flow Solution Procedures 

In order to determine the base flows about which perturbations 
are to be considered, two solution procedures were employed. These 
are (1) the quasi three-dimensional Euler equation solver B2DATL 
developed in references 12 and 13 which'solves equations (1) sub- 
ject to the boundary conditions depicted in figure 3, and (2) the 
two-dimensional small-disturbance transonic equation solver TSFOIL 
developed in reference 18 which solves equation (33) subject to 
boundary conditions appropriate to an unstaggered nonlifting 
cascade. While the B2DATL code provides the ultimate general 
capability for determining base flow solutions for the present 
perturbation study, run time and core storage requirements of that 
procedure limited the variety of perturbation applications which 
could be investigated. Notwithstanding the limitations of zero 
stagger, camber, and lift for cascade applications of the TSFOIL 
code, that procedure provides fast and accurate base solutions 
throughout the entire transonic range, including supercritical, 
choked, and supersonic flows. Figures 6 through 8 demonstrate the 
validity of that procedure to provide reliable base flow cascade 
solutions. Figure 6 displays results for the surface pressures 
predicted by TSFOIL for purely subsonic flow past two unstaggered 
nonlifting cascades with different pitch-to-chord ratios H/C and 
having biconvex blade profiles for several different thickness 
ratios. The plot on the left displays results for H/C = 0.75. 

The dashed line indicates the result of linear theory for t = 0.06 
and provides some idea of the magnitude of the nonlinear effects 
present even at this low a Mach number. The analogous results in 
the plot on the right for H/C = 2.0 indicate less severe but 
still appreciable nonlinear effects for this less confined flow. 
Both of these results display the need for a nonlinear calculation. 
Figure 7 exhibits the location of the sonic line and shock waves 
during the onset of choked flow as predicted by TSFOIL for a 
cascade of unstaggered, nonlifting parabolic-arc blade profiles 
with T = 0.07 and H/C = 2.0. We note the small Mach number range 
between the onset of supercritical flow and choking, and the move- 
ment of the shock wave downstream after choking has occurred with 
decreasing downstream pressure. Figure 8 exhibits the pressure 
distribution on the cascade centerline between the two blades. 

These results verify those measured and predicted for choked flow 
by Collins and Krupp (ref. 19) . Finally, figure 9 displays the sen- 
sitivity of surface pressure distributions to thickness ratio for 
this cascade for M = 0.80 as the thickness ratio varies from 
T = 0.065 to 0.075, and indicates again the high sensitivity of 
confined transonic flows. 
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3 . RESULTS 


3 . 1 Linear Perturbation Equation Method 

Theoretical results obtained from the linear perturbation 
equation method for various parameter perturbations are provided 
in figures 10 through 18. Initially presented are results based 
on the small-disturbance transonic equation formulation given in 
equations (31) to (36) . Those results extend through the tran- 
sonic regime and include subcritical, supercritical, and super- 
sonic flow situations. Also presented are subcritical perturba- 
tion results obtained using Euler equation base solutions. 

In figure 10, comparisons are displayed for perturbation 
solutions associated with thickness ratio changes of two unstag- 
gered nonlifting cascades composed of biconvex profiles in a sub- 
critical flow with free-stream Mach number M^ = 0.60. The plot 
on the left exhibits the comparisons for a pitch-to-chord ratio 
H/C = 0.75. 

The solid lines indicate the results obtained by using the 
base solution procedure (TSFOIL) for the three thickness ratios 
shown. The two solid line results for thickness ratios x equal 
to 0.06 and 0.08 are meant to be compared with those indicated by 
the dashed lines. The dashed results were obtained using the 
perturbation method based on solving equation (36) with an SLOR 
relaxation algorithm (CASCDE) subject to the boundary conditions 
for a thickness ratio perturbation given in Appendix B. The base 
flow solution for thickness ratio equal to 0.07 was used to 
evaluate the coefficients in equation (36) . The results with 
At = ± 0.01 were then added to the x = 0.07 base flow. 

The good agreement between the perturbation and base flow 
results is impressive considering that at this spacing the flow 
is relatively sensitive to changes in geometry. Numerical experi- 
ments show that this flow chokes at M^ just under 0.7. Further- 
more, the change in thickness ratio considered here of 1 in 7 is 
14 percent, which is not a small perturbation. Analogous results 
are presented in the plot on the right for H/C = 2.0. In this 
case, where the flow is less confined, the perturbation results 
are essentially identical with those predicted by using the base 
flow method. 

In figure 11, subcritical results are given for perturbation 
solutions associated with a change in pitch-to-chord ratio H/C 
of an unstaggered nonlifting cascade. This perturbation is 
fundamentally different from that associated with a Ax pertur- 
bation in its dependence on boundary data. As indicated in 
Appendix B and repeated in the sketch below, the primary charac- 
teristic of this perturbation is that the boundary conditions 
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involve the second derivative of the base flow potential along 
the dividing streamline and blade surface. Because of this, this 
particular perturbation problem provides a good test of the sen- 
sitivity of the perturbation procedure to boundary data involving 
base flow quantities. 

In figure 11, the dashed lines indicate the perturbation pre- 
dicted results for AH/C = ± 0.10. These were obtained by solving 
the perturbation equation (36) subject to the boundary conditions 
indicated above, and then adding those results to the H/C = 0.75 
base flow solution. Comparisons with the corresponding nonlinear 
TSFOIL solutions for H/C = 0.65 and 0.85 indicate reasonably 
good agreement for the perturbation change to a higher spacing 
(H/C = 0.75 to 0.85), and somewhat less satisfactory agreement for 
the change to a closer spacing (H/C = 0.75 to 0.65). However, the 
nonlinear TSFOIL results indicate a negative peak pressure coeffi- 
cient change approximately 50 percent larger for the H/C = 0.75 
to 0.65 change than for the H/C = 0.75 to 0.85 change. This 
clearly demonstrates the strong nonlinearity of the flow as the 
spacing is decreased. Furthermore, the nonlinearity is such that 
a AH/C decrease of 0.10 from the base flow at H/C = 0.75 is 
far beyond a linear change, while a corresponding increase is 
significantly closer to linear. 


A satisfying feature of the accuracy of these results is 
that, for a AH/C perturbation, the perturbation pressure coeffi- 
cients are composed of two terms; that is. 


(x,0) 

^u 



<|)^(x,0) + 


F 



with an analogous result for Cp.(x,H). For the results plotted 
in figure 11, F^(x) = - 4t = -0.28. Consequently, the perturba- 

tion pressure coefficient is the difference of two large quan- 
tities, so that the accuracy in the comparisons with the nonlinear 
TSFOIL results is even more impressive. 


Finally, we note that preliminary calculations of this per- 
turbation problem indicated a strong sensitivity to the accuracy 
of the (pyy distribution employed in the boundary data. It was 
found that grid clustering in the y direction near the bound- 
aries was essential in order to obtain accurate perturbation 
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results. This, of course, necessitated a somewhat finer grid and 
more accurately defined base flow than would normally have been 
used. 


In figure 12, the final subcritical results based on TSFOIL 
base solutions are presented for a perturbation in oncoming Mach 
number. In this instance, the boundary conditions for' the per- 
turbation are all homogeneous, but the governing equation 
(Appendix A) contains forcing terms on the right-hand side 
involving base flow quantities. We note that the perturbation 
problems involving changes in strecim sheet thickness b and radial 
divergence R, which have not yet been considered, involve the 
same basic problem formulation as this AMj perturbation; that 
is, homogeneous boundary data together with a governing differ- 
ential equation having forcing terms on the right-hand side. 

With regard to the comparisons shown in figure 12, we note 
an extremely strong nonlinearity of the flow as the Mach number 
is increased. In particular, the nonlinear TSFOIL solutions show 
that the negative peak pressure coefficient change from M^ =0.60 
to 0.50 is approximately 20 percent less than the Mj = 0.60 to 
0.65 change, which indicates a 120 percent difference per unit 
M^ change between increasing and decreasing Mj . Consequently, 
while the linear perturbation prediction falls somewhat short in 
predicting the M^ decrease, the behavior for an M^ increase 
is so strongly nonlinear that the linear method does not apply 
over a AMj =0.05 range. 

In all of the subcritical perturbation results provided in 
figures 10 through 12, the stability and convergence properties 
of the CASCDE solutions were excellent. No problems were 
encountered and convergence was usually obtained within 200-300 
SLOR iterations on a fine (103x21) mesh. 

Results for perturbations of supercritical cascade flows 
with TSFOIL base solutions are shown in figures 13 and 14. In 
figure 13, comparisons between the perturbation and base flow 
methods are given for a thickness ratio perturbation of the 
strongly supercritical flow past an unstaggered nonlifting cascade 
of parabolic-arc blade profiles with H/C = 2.0, Mj = 0.80, and 
base thickness ratio t = 0.075. A perturbation solution for 
At = 0.002 has been added to the t = 0.075 base solution to 
produce the dashed line result. Comparison with the nonlinear 
TSFOIL result for t = 0.073 indicates that the perturbation 
prediction is essentially identical with the nonlinear result 
until the shock wave is approached, at which point a discrepancy 
appears. Beyond the shock, agreement again becomes good. The 
reason for the discrepancy near the shock is that the conserva- 
tion form of the differential equation used to calculate the 
perturbation solution (eq. (36) ) essentially satisfies the appro- 
priate perturbation jump condition, but at the old (i.e., base 
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flow) location rather than new location of the shock. If the 
shock were somehow constrained not to move due to the flow pertur- 
bation, the present formulation would be valid. Thus, in the 
present formulation the proper shock condition is satisfied but 
at the wrong location. To remedy this, the perturbation problem 
must be supplemented by an additional condition at the base flow 
shock points which account for shock displacement. A similar 
condition holds true at sonic points, which are also displaced. 
However, since the flow is continuous at those points, the 
supplementary conditions there, although necessary to impose, 
would not be expected to have a large effect on the perturbation 
solution. 

In figure 14, analogous results are presented for a super- 
critical AH/C perturbation. In this instance, a significant 
discrepancy exists in the perturbation solution in the vicinity 
of the shock wave, where the perturbation solution exhibits a 
large oscillatory behavior. The cause of this discrepancy is 
clearly evident in the plot on the bottom of the figure, which 
displays the boundary data; that is, <j)yy along the blade surface. 
Extreme irregularities in the form of strong oscillations, absent 
in the subcritical case, are now present and located in the vicin- 
ity of the shock wave. Application of the supplemented perturba- 
tion shock point condition, as discussed previously, would no 
doubt improve the above comparison. However, the irregular behav- 
ior of the boundary distribution of <f>yy shown above may in fact 
be proper in the vicinity of the shock since large accelerations 
would be expected to occur. Should this be the case, it would 
be difficult to achieve accuracy for the Ah/C perturbation with 
the linear perturbation equation method for both supercritical 
and supersonic flows. Nevertheless, in spite of these irregu- 
larities in boundary data no convergence or stability problems 
were encountered in the CASCDE solutions for these flows. 

In figures 15 and 16, results are presented for perturbations 
of supersonic cascade flows with TSFOIL base solutions. Figure 15 
displays comparisons for a thickness ratio perturbation of 
At = 0.005 of an unstaggered nonlifting parabolic-arc profile 
cascade with Mj = 1.15, H/C = 2.0, and t = 0.70. The effect of 
not taking proper account of shock displacement is clearly evident 
in this example. The displacement of the bow shock, which is 
generally more sensitive to movement than surface shocks, is not 
well accounted for within the present formulation and may be the 
cause of the discrepancies on the blade surface. Corresponding 
results shown in figure 16 for an upstream Mach number perturba- 
tion, however, display quite good agreement. Further testing 
would be required to determine whether this is fortuitous or 
whether the Mj perturbation actually does become less sensitive 
and amenable to linear analysis at supersonic Mach numbers. 
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The final comparisons presented using the linear perturbation 
equation method are provided in figures 17 and 18. In these 
results, the base flow solutions are subcritical Euler equation 
solutions obtained from the B2DATL code (refs. 12,13). The cascade 
geometry used was identical to that employed in the TSFOIL solu- 
tions; that is, an unstaggered cascade of parabolic-arc blade 
profiles. In addition, in order to facilitate comparison with 
previous TSFOIL results, a small-disturbance perturbation formula- 
tion was used. Figure 17 displays comparisons for a thickness 
ratio perturbation of a nonlifting cascade with Mj = 0.6, 

H/C = 2.0, and t = 0.07. We note that in spite of the irregularity 
in the base flow solution near the leading edge, which is appar- 
ently due to the particular treatment of the sharp leading edge, 
and also the erroneous asymmetry of the solution about the blade 
midpoint, the perturbation procedure with At = ± 0.01 predicts 
excellent agreement with the base flow results at t = 0.06 and 
0.08. A similar comparison is given in figure 18, which displays 
results for a At = -0.01 perturbation about the same base flow 
conditions and geometry but now including lift, with the inlet 
flow angle 3i = 4°. Again, the agreement between base flow and 
perturbation predictions is excellent. 


3.2 Nonlinear Differencing Perturbation Method 

The initial impetus for application of the nonlinear differ- 
encing perturbation method was twofold: (1) to attempt to correct 

the previous deficiencies associated with shock displacement in 
the supercritical and supersonic perturbation problems discussed 
above, and (2) to determine whether this method would provide more 
satisfactory results for these more highly nonlinear flows than 
the linear perturbation equation method. Both of these charac- 
teristics are evident in the example selected for preliminary 
calculation by this method. This is the supercritical flow ini- 
tially presented in figure 9 for flow past an unstaggered non- 
lifting cascade of parabolic-arc blades of various thickness 
ratios with H/C = 2.0 and Mj = 0.80. The extreme sensitivity of 
the solution to small changes in thickness ratio is evident. The 
corresponding At perturbation solution via the linear perturba- 
tion equation method was provided in figure 13 and displays the 
typical shock displacement discrepancy. 

The source of difficulty as well as the remedy associated 
with defining a proper perturbation solution when shock waves are 
present and displaced by a perturbation is provided graphically 
by the sketches in figure 19. The sketch on the left illustrates 
as the shaded area the perturbation between two nonlinear super- 
critical solutions. We note that the perturbation is well defined 
and small everywhere except in the vicinity of the shock waves. 

In that region, due to shock movement, the perturbation becomes 
of the same order as the base flow solution and the entire 
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perturbation concept breaks down. The correction to this is 
indicated in the figure on the right in which coordinate straining 
is used to move the shock of the perturbed flow back to the base 
flow location. The result is that the perturbation now remains 
small and well defined everywhere. 

Initial results obtained by using such a coordinate straining 
technique have been recently reported (ref. 20) for small-dis- 
turbance flows past an isolated, nonlifting symmetric airfoil in 
which the shock wave that appears on the airfoil surface is 
assumed to be normal to the flow. We have extended that result 
to the cascade flows shown in figure 9. The results obtained by 
using coordinate straining together with differencing of two non- 
linear solutions are shown in figure 20. Here, the base flow 
solutions shown in figure 9 for x = 0.075 and 0.073 were dif- 
ferenced, after shock straining was employed, to define the unit 
perturbation solutions. Then, that perturbation solution was 
employed together with the base solution of x = 0.073 to pre- 
dict the solutions for x = 0.072, 0.070, 0.068, and 0.065. In 
each of the cases shown, the perturbation results (O) are meant 

to be compared with the dashed ( ) solutions which represent 

the nonlinear TSFOIL solution at the new thickness ratio. We 
note that in every case the agreement is outstanding, particularly 
in view of the large extrapolations involved. The only exceptions 
are several of the flagged (O') points designated on the figure as 
shock-capture points. Those discrepancies are not the fault of 
the straining procedure, but rather arise because of the shock 
smearing characteristic of the shock-capture procedure employed 
in TSFOIL. The coordinate straining procedure is derived on the 
basis of a sharp shock discontinuity. When shock-capture methods 
are employed for determining base flow solutions, the shock points 
will be somewhat smeared (refs. 18 and 21) as shown in figure 20. 
Corrections to this can be made directly by modifying the pertur- 
bation procedure in the vicinity of the shock either by smoothing 
or weighting the shock-capture points. 
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4. 


CONCLUSIONS AND RECOMMENDATIONS 


Theoretical analysis and associated development of computa- 
tional programs have been carried out in a preliminary application 
of perturbations methods to transonic flows in turbomachines. The 
theoretical analysis involved development of perturbation methods 
for determining first-order changes in the flow field due to 
variations in one or more geometrical or flow parameters for 
transonic flows on the blade-to-blade surface of a single blade- 
row compressor. The perturbation formulation is based on the 
quasi three-dimensional full potential approximation to the flow 
and is developed in a curvilinear coordinate system fixed to the 
blades. Included in the formulation are variations in stream 
sheet thickness and radius. Applications and results are pre- 
sented for selected two-dimensional cascade flows. 

Two different perturbation approaches were identified and 
studied. They are: 

• Linear perturbation equation method 

• Nonlinear differencing perturbation method 

With regard to the linear perturbation equation method, boundary 
value problems associated with a variety of perturbations were 
developed based on the full potential formulation. These include 
variations in maximum blade thickness, blade spacing, inflow Mach 
number, inflow angle, maximiam blade camber, and stagger. Also 
developed were the special forms of these boundary value problems 
based on the transonic small-disturbance formulation for maximiam 
blade thickness, blade spacing, and inflow Mach number. 

Perturbation solutions obtained with the linear perturbation 
equation method generally indicate that good results can be antic- 
ipated for blade geometry perturbations, such as blade thickness 
and angle of attack, for subcritical flows. This was found to be 
true both for perturbations of transonic small-disturbance flows 
as well as of the Euler equation flows developed in reference 12. 
Reasonable but less satisfactory results were obtained for per- 
turbation changes in overall quantities, such as blade spacing 
and free-stream Mach number, particularly when H/C was less 
than 1 and the flow was supercritical. Because these perturba- 
tions alter the basic character of the flow more rapidly, it is 
anticipated that they will prove to be more difficult to predict 
satisfactorily. Those results serve to point out the primary 
limitation of the linear perturbation equation method; that is, 
for near-critical transonic flows and blade spacings typical of 
modern compressor blade rows, the basic linear variation assump- 
tion fundamental to the technique is severely restrictive. 


25 



For the perturbation problem associated with blade spacing 
changes, where the boundary conditions as well as the differential 
equation involve base flow quantities, a sensitivity to base solu- 
tion accuracy was found. A direct remedy to this is to recast the 
perturbation problem in coordinates strained locally in such a 
fashion as to eliminate the base flow quantities from the bound- 
ary conditions. Although this procedure would complicate the 
perturbation differential equation somewhat, the sensitivity of 
the solution would be lessened and the overall perturbation solu- 
tion enhanced. 

Additionally, it was found that for flows with shock waves, 
a supplemental perturbation boundary condition is required at the 
shock wave. This is needed to account for shock displacement. 
Results obtained for supercritical and supersonic flows without 
this boundary condition are clearly unsatisfactory. 

The stability and convergence properties, however, of the 
linear perturbation equation code (CASCDE) that was developed to 
obtain finite-difference solutions of the various perturbation 
problems were very good for all the perturbations examined. This 
included the subcritical, supercritical, and supersonic cases 
considered, and held true in spite of strong irregularities in 
the base flow solution, such as displayed in the boundary data 
for the supercritical spacing perturbation (fig. 14) and the 
strong leading-edge inaccuracies in the Euler equation solutions 
for thickness (fig. 17) and lift (fig. 18) . Consequently, rapid 
and accurate converged solutions are definitely achievable with 
this method. 

In order to examine perturbation results which account for 
shock displacement, preliminary application of the nonlinear 
differencing method with coordinate straining was made to a super- 
critical cascade flow. The results obtained display outstanding 
agreement with corresponding nonlinear calculations . Aside from 
some slight discrepancies associated not with the perturbation 
technique or the straining but rather with the shock smearing 
character of the base flow procedure used (TSFOIL) , the method 
was able to predict from two strongly supercritical solutions 
satisfactory results all the way down to the limiting subcritical/ 
supercritical flow. 

On the basis of this preliminary study, it is evident that 
the two perturbation methods examined possess definite merit for 
application to certain transonic flow problems in turbomachines. 
Results obtained for selected perturbation problems indicate that 
these methods are capable of providing neighboring solutions of 
good accuracy on an economical basis. Further development is 
needed, however, to provide a computational tool of practical 
utility for a wide range of perturbations. Because each of the 
perturbation methods studied possesses distinct advantages, we 
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suggest that development of both is warranted. With regard to the 
linear perturbation equation method, inclusion of the shock dis- 
placement condition is needed and can be accomplished via coordi- 
nate straining. Similarly, for those perturbation problems whose 
boundary conditions involve base flow quantities, local coordinate 
straining should be used to eliminate that dependence and corre- 
sponding solution sensitivity. Finally, a variety of cases should 
be studied at flow conditions throughout the transonic regime to 
establish solution accuracy and range of validity. With regard 
to the nonlinear differencing method, in view of the success of 
this method in the initial application to a supercritical cascade 
flow, it appears to be an ideal technique for perturbing certain 
highly sensitive transonic flows. Consequently, it should be 
applied to a wide variety of perturbation problems and, in par- 
ticular, to strongly supercritical cases. In addition, two- 
dimensional straining should be incorporated in order to enable 
prediction of flow field as well as blade surface properties. 



APPENDIX A 


DERIVATION OF PERTURBATION DISTURBANCE POTENTIAL EQUATIONS 

This appendix outlines the deriv^ation of and provides expres- 
sions for the coefficients [A,B,C,D,E,F1 associated with the dif- 
ferential eqi^ation (11) governing the perturbation disturbance 
potential <J)^. The basis of the analysis is the nonrotating (fi=0) 
form of equation (6) given by 



a^) $, 


mm 



+ 



^m'^9^m6 


^ '^m d (Rb) 
Rb dm 


(A-1) 


where 


a2 = a2 + i(Y - 1) [U^ + - ($g/R)^l (A-2) 

and the subscript i denotes uniform inlet conditions. 

It is convenient to employ a disturbance potential <fi defined 
according to 


where 
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The perturbation coefficients [A,B,C,D,E,F] are determined 
by first expanding the potential 4> and the various parameters 
in the differential equations as follows: 
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By inserting equations (A- 3) through (A-9) into equations (A-1) 
and (a- 2) ,expanding the result/ neglecting terms of higher than 

first order in <|)\ Ag , AM , and A-^^^~, and then ordering the 

' / ^ 0 

resulting equations according to base flow ( ) and perturbation 

(<j)\ A8 , AMj, quantities, the following differential 
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APPENDIX B 

BOUNDARY CONDITIONS FOR PERTURBATION DISTURBANCE POTENTIALS 

This appendix contains the perturbation boundary conditions 
associated with various parameter perturbations of the quasi 
three-dimensional full potential base flow formulation described 
in Section 2.4.2. These include variations in maximum blade 
thickness, blade spacing, inflow Mach number, inflow angle, 
maximum blade camber, and blade stagger. In addition, special 
forms of these for maximum blade thickness, blade spacing, and 
inflow Mach number perturbations based on the transonic small- 
disturbance formulation for an unstaggered cascade are also 
presented. 

The coordinate system employed for the disturbance full 
potential formulation is sketched below. 


<i>T=0>B 





In this figure, F^(x) , Fj^(x) denote the dimensionless 
ordinates of the blade upper and lower surface and are given in 
normalized form by 


F^(x) = xf (x) + hg(x) (B-1) 

F^(x) = - xf(x) + hg(x) (B-2) 


where _x and_h represent the thickness ratio and maximum camber, 
and (f(x), g(x)) are the corresponding normalized thickness and 
camber distributions. The normalized quantities (x^,y^) are 
related to the blade spacing t and stagger angle X by 


t . , 

^ sin A 

y^ = ^ cos X 


(B-3) 
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and are the uniform inflow and outflow angles, respec- 

tively, of the absolute flow. It is convenient to define the 
quantities : 


= X + -2- 


^2^ ^ - T 


yi = T- 


(B-4) 
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X, 


0 -t 

x^ = X + etc. 

u The appropriate boundary conditions for both the base flow 
<p and the perturbation component for various perturbations 

are as follows: 

1. Blade thickness ratio perturbation 
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Periodic condition: 

- i j < X < 0 

0 

0 , 

<t> (x^, y^) = 

<t> (Xj, yj) 

<P^ (x^ , y^) = 

<P^ (Xj/ y^) 
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Cyclic condition : 

(Xj , Y^) = 

Yj) = 


1 < X < 1 + 5,. 


(j)°(x^, Ya) + 

y,) + 


(B-8) 


2 . Blade spacing perturbation 


t _ t At 

C C C 


Blade surface condition; 0 ^ x <. 1 
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Periodic condition ; - < x < 0 
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Continued on next page 


34 



(B-12) 
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3. Inflow Mach niimber perturbation 

Mj = M° + AM (B-14) 

Blade surface condition: 0 :< x ^ 1 


4)“(x^,Y,) 

= F^(x)cos 


- sin 3^ 


4>y(x^ry^) 

1 

= p (x) cos 
u 


- sin 3j 


4>y(x^,Y,) 

= 0 




4>y(X2,Y2) 

= 0 
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Periodic condition : - £ ^ < x < 0 

<l>° (Xj fYi) = <1>” (^2 fYg) 

<l>^(Xi,y^) = (B-16) 

Cyclic condition ; 1 < x < 1 + 

^'’(Xj/Yj) = (^'’Cx^/Yg) + 

4>^(Xj,Yi) = tl^Nx^/Yz) + (B-17) 
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4 . Inflow angle perturbation 


= e° + (B-18) 

Blade surface condition; 0 ^ x :< 1 


<}>y(x^,y^) 

1 0 

= F^(x) cos 3j “ sin 

< 


<^;(x2,y,) 
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Periodic condition; - < x < 0 

<!>'’ (Xj ,Yj) = (j)° 

= (J)^(x2,y2) (B-20) 

Cyclic condition ; 1 < x < 1 + 

<j)“(x^,yi) = 4)°(x2,y^) + A(j)“g 

4>^(Xj,y^) = ((^^(x^/y^) + (B-21) 

5 . Blade maximum camber perturbation 

h = h° + Ah (B-22) 

Blade surface condition ; 0 ^ x ^ 1 

(|)y(Xi,yi) = [-Tf ' (x) + h°g'(x)]cos 3j - sin 3^ 

= [Tf ' (x) + h°g'(x)]cos 3j - sin 3^ 

Continued on next page 


36 



<j)y(x^,yi) = Ahg' (x) cos 3^ 


= Ahg' (x) cos 3j 


Periodic condition ; - < x < 0 

4>° ,yj) = (()“ (x^ /y^) 

<P^ (Xj ,y^) = (x^/y^) 

Cyclic condition ; 1 < x < 

4»°(Xi,yj) = (j)“(x^,y^) + Acf)®^ 

4>^(Xj,yj) = <f)^(x2,y2) + 


(B-23) 
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6. Blade stagger angle perturbation 

X = X° + AX 

Blade surface condition ; 0 ^ x \< 1 
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Periodic condition ; - < x < 0 

<^°(xj,y;) = /(x",y°) 

cj)'(xj,yj) = 4>^(x°,y°) - [4>x^xJ,yJ) cos x' 


,0,0 0 . . , 0 

- 4>y(x2fy2) sin X 
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Cyclic condition ; 1 < x < 1 + 
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For an unstaggered cascade, the special forms of the boundary 
conditions for the perturbations given above follow directly. 
However, for an unstaggered nonlifting cascade, the additional 
symmetry of the flow provides a convenient alternative specifica- 
tion to the Dirichlet-type periodic and cyclic boundary conditions 
previously derived. This is, of course, based on the fact that 
the horizontal boundaries extending fore and aft of the blades, 
due to the flow symmetry, become dividing streamlines so that a 
Neumann condition of no normal flow can be specified. The basic 
form of the boundary conditions for the disturbance potential <p 
then become 
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This leads to the following boundary conditions for 
perturbation problems associated with thickness ratio x 
ratio H/C, and free stream Mach number changes. 

1. Thickness ratio perturbation 

0 . 

X = T + Ax 


Blade surface condition: 0 ^ x :< 1 


1 




4> 


I 

y 




Axf ' (x) 
Axf ' (x) 


Dividing streamlines ; -SL^ < x < 0, 


1 < X < 1 + il 
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2. Spacing ratio perturbation 
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(B-36) 


3. Inflow Mach nximber perturbation 

M^ = M° + AMj 


Blade surface condition: 0 ^ x ^ 1 
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Dividing streamlines : 


<x<0, I<x<l+Jl 2 
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For the Mach number perturbation, the perturbation boundary 
conditions are homogeneous but the differential equation is not. 
For the transonic small-disturbance approximation, instead of 
the perturbation equation given by equation (36), we have instead 
the following equation. 



m"(Y+1)4> 
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APPENDIX C 


LIST OF SYMBOLS 

a local speed of sound/ m/sec; eq. (7) 

A coefficient of 4>xx perturbation disturbance 

potential differential equation (11) ; defined by eq. 
(A-13) 

b stream sheet thickness, m 

B coefficient of (p^y in perturbation disturbance 

potential differential equation (11) ; defined by eq. 

(A- 14) 

C blade chord, m 

C coefficient of <|)yy in perturbation disturbance 

potential differential equation (11) ; defined by eq. 
(A-15) 

C pressure coefficient, (p-p )/4 P 

1 « 1 1 

D coefficient of 4>^ in perturbation disturbance 

potential differential equation (11) ; defined by eq. 
(A-16) 

E total energy per unit mass, m^/sec^; eq. (2) 

relative total energy per Unit mass, m^/sec^; eq. (2) 

E coefficient of (j)y in perturbation disturbance 

potential differential equation (11) ; defined by eq. 
(A-17) 

r nondimen sional ordinates of blade thickness; normalized 

by tC 

nondimensional ordinates of the blade upper and lower 
surfaces, respectively; normalized by C 

F right-hand-side of perturbation disturbance potential 

differential equation (11) ; defined by eq. (A-18) 

g normalized ordinates of blade mean camber line; 

noarmalized by hC 

h blade maximum camber ratio; i.e., maximum height of 

mean camber line normal to chord divided by chord 
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total enthalpy per imit mass, m^/sec^; eq. (2): also, 

blade spacing for nonstaggered cascades , m 

relative total enthalpy (rothalpy) per unit mass, 
m^/sec^; eq. (2) 

transonic small-disturbance similarity parameter; 
eq. (34) 

nondimens ional length denoting distance from blade 
leading edge to upstream boundary, normalized by C 

nondimensional length denoting distance from blade 
trailing edge to downstream boundary, normalized by C 

curvilinear blade- fixed coordinate denoting downstream 
direction, m 

absolute inlet Mach ntimber 
static pressure, Newton/m^ 
uniform inlet velocity, m/sec 
uniform outlet velocity, m/sec 

approximate perturbed flow solution for arbitrary flow 
quantity, eq. (3) 

base flow solution for arbitrary flow quantity, eq. (3) 

linearized perturbation solution per unit perturbation 
change for arbitrary flow quantity, eq. (3) 

stream sheet radius of blade-to-blade surface, m 

blade spacing in 0 direction, m 

relative velocity component in blade-fixed m-coordinate 
direction , m/sec 

relative velocity component in blade-fixed 0-coordinate 
direction , m/sec 

relative velocity-component normal to blade surface in 
blade-fixed coordinate system, m/sec 

nondimensional blade- fixed coordinates related to the 
(m,0) coordinate system by eq. (8) , normalized by C 



nondimensional blade spacing quantities defined by 
eq. (12) 

nondimensional transonically scaled y-coordinate , 
defined by eq. (31) 

absolute inflow angle, radians 

absolute outflow angle, radians 

ratio of specific heats, equal to 1.4 for air 

amplitude of perturbation of geometric or flow parameter 

absolute angular coordinate, defined in figure 2; 
radians 

blade setting or stagger angle, radians 
density, kg/m^ 

blade thickness ratio; i.e., maximum blade thickness 
normal to chord divided by chord 

dimensionless disturbance velocity potential, normalized 
by Cq^; eq. (8) 

dimensionless base flow disturbance velocity potential, 
normalized by Cq^ ; eq. (10) 

dimensionless perturbation disturbance velocity poten- 
tial, normalized by Cq^ ; eq. (10) 

dimensionless transonic small-disturbance velocity 
potential, normalized by Cqj ; eq. (32) 

dimensionless base flow transonic small-disturbance 
velocity potential, normalized by Cq^ ; eq. (35) 

dimensionless perturbation transonic small-disturbance 
velocity potential, normalized by Cq^ ; eq. (35) 

jump in the disturbance velocity potential at the blade 
trailing edge 

dimensional velocity potential for the complete flow 
defined in the absolute (nonrotating) reference system; 
eg* (5) , m/sec 

rotational speed, radians/sec 



Subscripts 

u refers to blade upper surface 

% refers to blade lower surface 

^ denotes uniform inlet conditions; also denotes upper 

boundary of solution domain 

^ denotes uniform outlet conditions; also denotes lower 

boundary of solution domain 

Superscripts 

“ denotes base flow quantities 

^ denotes perturbation quantities 
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Figure 3.- Boundary conditions for the quasi three-dimensional Euler equation 
formulation on the curvilinear blade-to-blade stream surface for a 

system rotating with the blades. 



Figure 4.- Boundary conditions for the quasi three-dimensional disturbance 
full potential formulation on the curvilinear blade-to-blade stream 
surface for a coordinate system rotating with the blades. 




where [G,H, I, J,K,L] are in general functions of '^yy ^xx^ 


Figure 5.- General form of the boundary conditions 
for the perturbation disturbance potential, c{i\ 






(a) Note small Mach number range between onset bf transonic flow and choking. 

(b) Choked conditions represent a singular point for the perturbation 

procedure - can only perturb conditions downstream of the limiting 
characteristic . 
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Figure 7.- Location of sonic lines and shock waves during onset of choked 
flow as predicted by TSFOIL for an unstaggered, nonlifting cascade 
of parabolic-arc blade profiles with T =0.07 and H/C = 2.0. 






Figure 9.- Surface pressure distributions on an unstaggered nonlifting 
cascade of parabolic-arc blade profiles for various thickness 
ratios with M, =0.80 and H/C =2.0 by TSFOIL. 




Parabolic- Arc Blades t = 0.07, Unstaggered Nonlifting Cascade 

Base H/C = 0.75, Perturbation AH/C = ±0.10 


IJ /O — A ^ C 
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Figure 11.- Subsonic H/C perturbation solutions obtained 
using TSFOIL base solution and CASCDE perturbation code. 



Figure 12.- Subsonic Mj perturbation solutions obtained 
using TSFOIL base solution and CASCDE perturbation code. 






Figure 13.- Supercritical perturbation solution obtained 

using TSFOIL base solution and CASCDE perturbation code. 
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Parabolic- Arc Blades x = .075, Unstaggered Nonlifting Cascade, M = 

Base H/C = 2.0, Perturbation AH/C = 0.2 



Figure 14.- Supercritical H/C perturbation solution obtained 
using TSFOIL base solution and CASCDE perturbation code. 
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Figure 15.- Supersonic 'T perturbation solution obtained 
using TSFOIL base solution and CASCDE perturbation code. 
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Figure 16.- Supersonic perturbation solution obtained 

using TSFOIL base solution and CASCDE perturbation code. 
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Figure 17.- Subsonic t perturbation solutions obtained 
using B2DATL base solution and CASCDE perturbation code. 




Figure 18, - Subsonic perturbation solutions obtained 

using B2DATL base solution and CASCDE perturbation code. 



Ordinary perturbation solution - 
Obtained by differencing two nonlinear 
solutions given in physical coordinates 

X. 



Shaded area represents 
perturbation solution. 

Note large perturbations 
resulting from shock 
movement . 


New perturbation solution - 
Obtained by differencing two nonlinear 
solutions after introducing strained 
coordinates x ' such that two shock 
locations (and leading and trailing 
edges) are brought together. 



Shaded area represents 
perturbation solution. 

Note large perturbations 
associated with shock 
movement have been removed. 
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Figure 19.- Perturbation shock displacement and proper definition 
of perturbation solution via method of strained coordinates. 
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